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We investigate the strong deflection limit of gravitational lensing by a Schwarzschild black hole 
embedded in an external gravitational field. The study of this model, analogous to the Chang & 
Refsdal lens in the weak deflection limit, is important to evaluate the gravitational perturbations on 
the relativistic images that appear in proximity of supermassive black holes hosted in galactic centers. 
By a simple dimensional argument, we prove that the tidal effect on the light ray propagation mainly 
occurs in the weak field region far away from the black hole and that the external perturbation can 
be treated as a weak field quadrupole term. We provide a description of relativistic critical curves 
and caustics and discuss the inversion of the lens mapping. Relativistic caustics are shifted and 
acquire a finite diamond shape. Sources inside the caustics produce four sequences of relativistic 
images. On the other hand, retro-lensing caustics are only shifted while remaining point-like to the 
s_i ■ lowest order. 
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I. INTRODUCTION 



Gravitational lensing in the Strong Deflection Limit (SDL) is emerging as a promising tool for black hole investi- 
gation. It is well known that photons making one or more complete loops around the black hole before approaching 
the observer produce two infinite series of images very close to the shadow of the black hole 0,0, 0- However, this 
(T) • fact became an astrophysical breakthrough only when the super-massive black hole supposed to be hosted in the 
radio-source Sgr A* in the Galactic center was proposed as an ideal candidate for detecting relativistic images of 
sources passing behind it . SDL investigations strongly benefitted from expansion techniques of the deflection angle 
[j| , which were then generalized to arbitrary spherically symmetric metrics 01. Theoretical prop erties of SDL lensing 
are now well assessed for a large class of black holes, since effects of charge U0, rotation [i^.llO| and several extended 
O" 1 theories of gravitation HIHElEQIlEIISIHEi have been considered - 

5^ \ On the observational side, lensing effects are expected to play a significant role in all high-resolution observations. 
bS)- However, a clean detection of a deeply relativistic image of some source is by no way an easy task, since Sgr A* must 
be resolved at least at microarcsecond level, and a source is needed with a surface brightness much higher than Sgr A* 
itself in the spectral band in which observations are lead. The development of Very Long Baseline Interferometry 
(VLBI) techniques in different bands of the spectrum is fast approaching the microarcsecond resolution required to 
distinguish such relativistic images (for an updated discussion on observational perspectives of relativistic images, 
see e.g. Ref. 0). As discussed in Ref. 0, the best prospect comes from relativistic images of low mass X-ray 
binaries around Sgr A* 20], which should be detected by future X-ray interferometry missions such as MAXIM 
(http:/ /maxim. gsfc.nasa.gov). 

Supermassive black holes should be hosted in galactic bulges and the surrounding environment could perturb the 
space-time of the main lens. As the black hole is part of a larger system, the gravitational field around it can not be 
exactly spherically symmetric. In some cases, large clusters of stars or secondary black holes are known to lie close to 
supermassive black holes. In this paper, we want to discuss the effect of external perturbations on the trajectories of 
light rays that are highly deflected. In the weak deflection limit, lenses with perturbed symmetry are usually modelled 
as quadrupole lenses |2lJ . If the gravitational field of the perturber changes very little over the relevant length scale of 
the black hole, it is enough to expand the potential of the external macro-lens up to the first non trivial term, i.e. the 
quadratic term. The perturbed Schwarzschild lens in the weak deflection limit is known as Chang-Refsdal lens after 
its first investigators j2^,|2^|. It was also recently revisited in Ref. j^. Here, we offer a complementary treatment in 
the SDL, which may be relevant to clarify the role of external fields in the formation of relativistic images. The paper 
is structured as follows. In Section \n\ we introduce the general lens equation for an unperturbed Schwarzschild black 
hole and specify it for the two cases of standard lensing and retro-lensing. In Section ITTT1 we analyze the perturbation 
on the deflection angle exerted by an external weak field, showing that it mainly occurs in the weak field regime. 
In Section IIVI we analyze the perturbed lens equation in the standard lensing geometry and in the retro-lensing 
configuration, finding critical curves, caustics and images. Section ]V\ is devoted to some final considerations. 
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FIG. 1: Generic lensing configuration. 



II. THE SCHWARZSCHILD LENS IN THE STRONG DEFLECTION LIMIT 



If the lens is spherically symmetric, as in the case of a Schwarzschild black hole, the motion of the photon takes 
place on a single plane. Then the lens equation relating the angular position of the source f3 and the angular position 
of a lensed image 9, relative to the position of the lens in the sky, can be found by simple trigonometry, looking at 
Fig. [U 

We denote by Dl, Dg and -Dls the distance between observer and lens, observer and source, and lens and source 
respectively. They are related by Carnot theorem D 



LS 



Z?g — 2DlD$ cos fl. We also introduce the angle 7 as 
the angular position of the source as seen from the lens. The angle 9 is the emission angle of the light ray from the 
source w.r.t the lens. The deflection angle q;bh is the angle between the asymptotic incident and outgoing directions 
of the photon. 

The impact parameter u can be taken as the distance between the lens and the incident direction or the outgoing 
direction. The two distances are indeed the same by time-reversal symmetry. This implies the equality 



Dl sin 9 = Dls sin 9. 

Let us consider the triangle OLS. Indicating the angle OSL by the sum of the internal angles gives 

C + (3 + (7T - 7) = 7T. 

Considering the triangle OCS, the sum of the internal angles gives 

(0-/3) + (3-O + (7t-obh) =*•■ 
Summing the two equations, we get the relation 

7 = 9 + 9 - qbhi 



Ds sin (3 = Z?ls sin 9 cos(dBH — 9) — Dls cos 9 sin(dBH 
Using Eq. (JJJ to eliminate 9, we finally get the general form of the lens equation 



Ds sin (3 — sin 9 cos(q:bh 



D£ s — Dl sin 2 9 sin(o!BH 



(5) 



(6) 



This equation is more general than the lens equation derived in Ref. Q , which was derived assuming that C belongs 
to the lens plane. That assumption is only valid for small «bh, making that equation inadequate for retro- lensing or 
other cases when the incident direction differs very much from the outgoing direction. Our lens equation is completely 
general and can be used to describe any configuration. 
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If the lens is an unperturbed Schwarzschild black hole of mass m, the deflection angle takes the classical weak 
field form Am/u for large impact parameters u, where we are using geometrized units (gravitational constant G = 1, 
speed of light in vacuum c = 1). For small impact parameters, the deflection angle grows more and more, diverging 
at u — u m — 3\/3to- The calculation of the exact deflection angle as a function of u has been first carried out by 
Darwin analytically in terms of elliptic integrals. A good approximation, valid for small impact parameters (large 
deflection angles), is the so-called Strong Deflection Limit (SDL), firstly introduced by Darwin himself and lately 
generalized by Bozza to any spherically symmetric black hole . It reads 

d B H ^ -a In I 1 ] + b, (7) 

For a Schwarzschild black hole, a = 1, b = —0.40002. Note that 6* m ~ u m /D^ corresponds to the angular radius 
of the apparent shadow of the black hole, which is typically very small in all cases of interest. For example, it is 
expected to be ~ 23 //as for the supermassive black hole at the center of our Galaxy. This formula is strictly valid 
for DlSjDl ^ 2m. For sources or observers very close to the black hole it needs to be corrected. 

As shown in several papers [ll|3j,|4j,|5|,|3,|23,|27j, relativistic images formed by light rays experiencing large deflections 
are prominent only when the source is in proximity of a caustic point. The caustic points are behind the lens and 
in front of the lens. The first class of caustic points is relevant for standard gravitational lensing configurations, in 
which the source is almost aligned behind the lens. The second class is relevant for retro-lensing configurations, where 
the source is on the same side of the observer and the light rays must turn around the black hole in order to hit the 
observer. In the following, we shall specify the general lens equation for sources close to these two classes of caustic 
points. 



A. Lens equation for standard lensing geometry 

Considering sources very close to the optical axis behind the black hole, we have /?, 7 <C 1, while the relation 
L\ S> 2m implies 9,9 <§; 1. Then, Eq. |@}, interpreted as an equation valid modulo 2tt, implies &bu = 2nir + Ad™, 
with Aa^ H « 1. In consequence of this, all trigonometric functions in Eq. © can be expanded for small angles. 
Taking the first order terms of these expansions, we get 

^9-^Aa™, (8) 

where we have also used that Dg ~ L\ + Dls- This is the lens equation for standard lensing configurations. The weak 
field lens equation is included for n = 0, i.e. when the deflection angle is very close to zero. Then Ad™ coincides 
with q;bh- On the other hand, for small impact parameters, the deflection angle grows very large, making the photon 
complete one or more loops around the black hole. In such cases, in the lens equation one has to consider the offset 
Ad^ H = dBH — 2mr. 

If the source is behind the lens, we have images for each integer n, corresponding to a deflection angle that is close 
to a multiple of 2tt. In consequence of this, we are mainly interested in those values of 9 such that cSbh is a multiple 
of 2ir. Indicating these values by 0°, we have a_ei/(L ) L^) = 2mr. From Eq. (7J, we get a simple expression for #° 

9° n ~6 m (l + e n ), e n = e^ 2n ^/\ (9) 
Then, letting e= (9 — 9^)/ 9^ and expanding the off-set around 6*°, we get 



The lens equation reduces to 



where 



/^C(l-i), (11) 



e n = ~T\ -V™- \ 12 1 



Dls a 

and we have already taken into account that e,tj«l to get rid of terms of order e 
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Images on the same side of the source solve Eq. 1)120 for a positive (3. The solution of the same equation with the 
source placed in —j3 gives the position of the images on the opposite side. If j3 — 0, the source is lensed in an infinite 
sequence of relativistic Einstein rings at 9 = 9 n (l + eJJ). For ^ 0, two infinite sequences of relativistic images form 
on opposite sides of the black hole at 

Note that for any relativistic image e » j^-9 m . Thus, unless Z?ls ~ 2m, we have e <C 1, so that 6*° already 

represents an excellent approximation to the position of the n-th relativistic image. 



B. Lens equation for retro-lensing geometry 



The treatment presented up to now only covers standard lensing geometry, with a source almost aligned behind 
the black hole. The modifications needed to account for the case in which the source is in front of the lens are very 
simple. In fact, redefining 6 n by the equation aBn(Di,6 n ) = (2n — 1)7T, we get 



_ e [6-(2n-l)ir]/o_ 



9° n =e m (l + e n ), 

The expansion of Eq. © for small angles and dsH = (2n — l)n + Ad™ leads to 



D s J D s 



Ad 



BH 



(14) 



(15) 



which differs from Eq. © for two signs coming from the expansions of the trigonometric functions and from the fact 
that now Ds = -Dl — £>ls- Since the magnification privileges sources close to the black hole, here we only consider 
Dls < Dl, while Ref. ^2^ considers the other case Dls > -Dl, with a slightly different expression for the lens equation. 
The expression for Ad™ remains the same given in Eq. 1 )10) 1 . so that the lens equation can be written as 



/3~- 1 



Ds J 



where now 



R _ 



Ds 

Dls a 



J 



(16) 



(17) 



The relativistic Einstein rings form at e = e^. For j3 ^ 0, the relativistic images form at 



i±A 

6° 



D s J 



(18) 



III. DEFLECTION ANGLE IN A PERTURBED SCHWARZSCHILD GEOMETRY 

We assume that the perturbation is weak, so that the gravitational potential tpQ induced by the tidal field is small 
throughout the photon trajectory. The order of the curvature radius TZ of the external field at the black hole position 
is then determined by the scale of the variations in the gravitational potential induced by the perturber. Indicating 
this scale by b, we can thus establish the relation 

Tl- 2 w dtdjipQ pb 6 -2 |pq|. (19) 

If the scale of the variations b is much larger than the Schwarzschild radius of the black hole 2m, this relation 
automatically implies 1Z 3> 2m. For example, if the external field is due to a second point-mass, we have <pq = 
— 2m2/r2, where is the distance of the generic point from the second mass. In this case, the scale of the variation of 
such a potential around the first black hole is dictated by the distance b between the two black holes and 1Z~ 2 w 7712/6 3 , 
coherently with what we know from the curvature invariants of the Schwarzschild metric. 

Coming back to the general situation, the metric of a spherically symmetric black hole distorted by an external 
gravitational field can be calculated with an expansion in inverse powers of the local radius of curvature 1Z of the 
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external spacetime. The first corrections, proportional to 1Z~ 2 , have been calculated by Alvi |28(, who found an 
approximate form of the metric of a binary system of black holes. In particular, he divided the spacetime in four 
regions: two of them surround the two black holes, an intermediate region is then encircled by the asymptotic region 
far from the two black holes. For each region, he found an approximated expression for the metric, nicely matching the 
metrics of the neighbor regions in the overlap zones. These approximate solutions were proposed as initial conditions 
for numerical programs designed to find the best approximation to the whole spacetime. Recently, Poisson [2j| gave 
the general expression for the tidal distortion of a Schwarzschild black hole embedded in an arbitrary spacetime up 
to terms proportional to 1Z~ 3 . In his work, the external spacetime determines the tidal fields through the electric and 
magnetic components of the Weyl tensor. These fields are then multiplied by radial functions to be determined solving 
the Einstein equations. The metric found by Poisson reproduces Alvi's metric when the tidal fields are generated by 
a second Schwarzschild black hole. What matters for us, as we shall see in a while, is that the first corrections close 
to the black hole are always of order of TZr 2 . 

In an asymptotically flat spacetime the deflection vector is simply given by the difference between the incoming 
and outgoing tangent vectors to the photon path. Integrating the geodesies equation, we get 

T\^x u dr, (20) 

where T % ^„ is the affine connection and r is an affine parameter. The integral can be split into three pieces: one 
covering the approach to the black hole, where the effects due to the black hole and the perturber are both weak; a 
second piece where the photon winds around the tidally distorted black hole under its strong gravitational field; a 
third piece where the photon goes away from the black hole and the deflections due to black hole and perturber are 
again weak. Denoting by — rstrong and Tstrong the values of the affine parameter where the transitions between the 
different regimes take place, we have 



a = d in + "Strong + "out, (21) 

^Strong 

j Y\ v x»x v dT, (22) 



-TStro 

din — 



— oo 

TStroi 



"Str 



ong 



"out = 

TSt 



J r\„x»x u dT, (23) 

"Strong 
OO 

J Y\ v x»x v dr. (24) 



The first and the third piece can be calculated in the weak field regime, where T l ^x fi x L/ ~ 2(Vj_<^)\ with W±ip 
being the projection of the gradient onto the plane orthogonal to the direction of the light ray and ip the Newtonian 
potential, given by the sum of the black hole and perturber contributions (p — </?bh + VQ- The central piece, dstrong, 
can be expanded in inverse powers of the tidal radius 1Z. Following Poisson |29j. the first order correction dg^ rong is 
of order of 1Z~ 2 . We can thus write 

— ^Strong OO 

2 J V±<pBHdT + c4°i ong + 2 J S7±ip B iidT 



a 



— OO 
-TStrc 



+ 2 J V±VQdT + a««g + 2 J Vxyqrfr. (25) 



rstrong 



The first three pieces reconstruct the unperturbed Schwarzschild deflection vector qbh, related to the scalar deflection 
angle introduced in the previous section by 



|&bh| =2 



. Abh — 2n7r 
sm 



(26) 



Note that the deflection vector, as defined by Eq. H2()(l , is independent of the number of loops performed by the photon 
around the black hole, since it only refers to the asymptotic directions. 
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Now let us consider the terms introduced by the perturber. The perturbation in the metric of a tidally distorted 
black hole is of the order of 7^ 2 , i.e. Sg^ w (r/TZ) 2 29]. Then, ST 1 ^ w {r/TZ 2 ). The contribution to the deflection 
angle due to the external perturbation is then given by an integral over the strong held region, whose size is governed 
by the Schwarzschild radius of the black hole 2m, of the perturbation to the afhne connection in the tidally distorted 
Schwarzschild metric, 



Strong 



2 m 



K 2 



dr 



2m 



(27) 



This is what expected on a dimensional analysis since o^trong ^ TZ 2 and, apart from the local curvature radius of 
the external field providing the scale relevant for tidal perturbations, the only other length is 2m itself. 

On the other hand, the two weak field contributions contain integrals of Vj_(£>q ~ b~ 1 ipQ and cover the whole 
photon trajectory, safe for a small piece whose length is of order 2m. Since the scale of variations b determines the 
size of the region where ifQ is relevant, the two weak field contributions are of order 



— ' strong oo 

2 J V±<p Q dT + 2 J X7 ± tp Q dT 



\<PQ\ 



(28) 



Recalling Eq. I|19[) . we see that the contribution to the deflection due to the perturber from the weak field region is 
much larger than the strong field contribution. 

We can thus conclude that the perturbation of the trajectory of a photon passing close to a black hole due to an 
external field is dominated by the weak field contribution and that the deviations arising in the strong regime can be 
safely neglected. If we replace the two integrals in the range ] — oo, — rg tr ong] U [jstrongj oo[ by the integral covering the 
whole dominion ] — oo, oo[, we commit a relative error of order 2m /b <C 1. We can finally express the whole deflection 
vector as 



a = ckbh + c*c 



oo 

oq ~ 2 J \7 ± ip Q dT. 



(29) 
(30) 



So we have recovered a superposition principle between the strong deflection due to the black hole and the weak 
deflection due to the external field. This has been possible since the perturbation of the photon trajectory arising in 
the strong field regime is very small compared to that arising in the weak field regime. As we shall see in the following 
sections, the consequence of this crucial fact is that the lens equation preserves a very simple structure, allowing a 
straightforward analysis of the effects of the external field on the relativistic images. 



IV. PERTURBED SCHWARZSCHILD LENS 



A. Standard configuration 



Since the strong deflection angle due to the black hole, c*bh, and the weak deflection angle due to the external field, 
ckq, can be added, the lens equation will take a very simple form. However, the lens equation (jSJ and its specifications 
in the standard JSJ) and retro-lensing l|15|) geometries can not be directly generalized to the presence of an external 
field. The superposition principle was derived assuming that total deflection occurs by steps. After the photons wind 
around the black hole, they quickly reach the region where the only weak gravitational field is effective and eventually 
approach the asymptotic trajectory to the observer with impact parameter u. In other words, the impact parameter 
mbh which enters in the expression for the deflection angle due to the black hole in Eq. Q differs from the impact 
parameter u as seen from the observer after the action of the perturber. Equation © was explicitly derived assuming 
«bh = Difi, which is no longer true when an external deflection is present, and can not be directly applied. 

Let us define a coordinate system centered on the black hole and with the z-axis oriented towards the observer. 
Then let us consider a photon coming out from the black hole with an impact parameter ubh- We shall indicate its 
projection on the lens plane by the vector «bh- As the photon travels towards the observer, it will still be deflected 
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by the external field. Its trajectory can be thus parameterized as 

TO TO 

x(t)=x +x (t-t )-2 J dr' J V^qdr", (31) 

T T' 

where xo = (0,0, Dl) is the observer position; xo = (— 0i, — 1) is the tangent vector when the photon reaches 
the observer (to first order in 6) and to — -Dl is the value of the affine parameter in which the photon reaches 
the observer. Note that here we are only including the weak field shift due to the external field, as the black hole 
deflection is nearly ineffective in the weak field region. We will come back later to this point for its full clarification. 
Furthermore, we are only interested in how the perturber affects the relation between the asymptotic trajectory and 
the impact parameter at the black hole position. 

The integral in Eq. I|31|) can be calculated along the unperturbed trajectory xS°>(t) = (u^ H , u BH , t). Moreover, 
the integrand is significant only for r < b so the result does not change if we push the upper integration limit to to 
infinity. So, evaluating x(r) at r = and identifying «bh with x(0), we finally have 

tt BH = x(0) = D h 6- Au , (32) 

OO OO 

Au = 2 / dz I S7±Lp Q (uf H ,u% H ,z")dz". (33) 



o 



Since u = D^O, we see how the presence of an external perturbation modifies the relation between angles and impact 
parameters. The lens and the relativistic images it produces are shifted by Auq/D^ in the observer sky w.r.t. their 
true position relative to the black hole. 

Let us estimate the order of magnitude of the shift term in Eq. (|33|l . By definition, Auo is the spatial shift due 
to the deflection action of the perturber between the black hole and the observer. As before, the transverse gradient 
introduces a factor b" 1 , while the two integrals multiply the result by b 2 . We thus have |Auo| ~ fqb. This quantity 
might be comparable with ubh itself. To understand this, as perturbing body let us consider a second black hole 
with mass mi at distance b. Then |Auo| ~ 2rri2 and if the second black hole has a mass similar to the first one, 
ubh ~ 2toi « 2m 2 ~ |Au |. 

The reader may wonder why in the standard weak field limit we do not need to evaluate this kind of shift of the 
lens plane. Indeed, there is a big difference. In fact, in the weak field, the scale of variations in the deflection angle is 
what we have called b and is very large if the photons pass very far from the black hole. Then, a shift of order (pqb in 
the lens plane is absolutely negligible on the scale of the deflection b. It is practically a higher order term in a weak 
deflection expansion. For relativistic images, where we need a resolution of the order of the Schwarzschild radius 2mi, 
this shift of the apparent lens position is of the same order of the relative positions of the images. 

Another point to clarify is the fact that we are assuming that the whole deflection by the black hole takes place in a 
region of order 2m\, neglecting any weak field deflection by the black hole. Indeed, the weak field deflection induced 
by the black hole is governed by the transverse gradient of its potential. Yet the photon is coming out with a very 
small impact parameter, so the transverse gradient is of order 2ttiiWbh/^ 3 ~ (2mi) 2 /& 3 while the transverse gradient 
of the potential of the perturber is </?q/&, so that the first is of higher order in the gravitational potential w.r.t. the 
second. So, we can safely assume that the photon coming out from the black hole reaches the asymptotic trajectory 
in few Schwarzschild radii and suffers no more bending by the black hole. This can be explicitly verified plotting the 
trajectory of a photon experiencing a very large deflection. 

To get the lens equation, we must relate 6 to (3. This can be done adding the information on the photon path 
between the source and the black hole in analogy to Eq. H32|) . In fact, the photon trajectory from the source to the 
black hole can be parameterized by 

_/ 

T T 

x(r) = x s + x s (r - r s ) - 2 J dr' J V^qdr", (34) 

Tg ts 

where xs = (Ds0x, Dgfa, — -Dls) is the source position; xs = (ai — 6\,cx2 — $2,1) is the tangent vector when the 
photon leaves the source (to first order in 6 and a) and rs ~ — Z?ls is the value of the affine parameter in which the 
photon leaves the source. Note that the deflection vector entering here is the total deflection given by Eq. 129|) . 

Let us determine the point where the photon hits the lens plane. Since we are neglecting the effects of the external 
field when the photon winds around the black hole, the motion around the black hole is still planar. Then, the 
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incidence point coincides with the exit point mbh- We thus set x(0) = «bh i n Eq. I|34|l . so that 

«bh = x(0) = D s (3 + D LS (a -9)- Au s (35) 



z 

Au s = 2 J dz 1 J V ^ip Q (u™,u™,z")dz". (36) 

— oo — oo 

Now, eliminating mbh between Eqs. <|32ll and (|35|l . we get the perturbed Schwarzschild lens equation in the standard 
geometry 

P = 8 - TT [«bh + olq] + . (37) 

Apart from the presence of the perturber deflection angle ckq, we have the two shifts Aug and Auo. By their 
definition, it is easy to verify that 

Aw s - Au = D LS a Q - 2 f dr' J V ± tp Q dr" . (38) 

The total shift is related to the discrepancy between the position of the center of the perturber and the plane of the 
black hole. It can be effective if the external field mostly intervenes before or after the black hole and it is of higher 
order and can be neglected in the lens equation only if the distance between black hole and perturber along the line 
of sight is of the same order of the Schwarzshild radius of the black hole. If the deflector potential has the symmetry 
<p Q (uf H , uf H , z) = ip Q (uf H , uf H , -z), then the shift is null. 

The perturbed lens equation l|37|l reproduces the vector generalization of Eq. (JSJ when the external field vanishes. 
The black hole deflection vector qibh m standard lensing reduces to the offset Ad n , introduced in Section [HI in its 
obvious vector generalization. Needless to say, both deflections Aa n and &q must be evaluated in Ubh = u Auq- 

At this point it is worth to eliminate 9 throughout the lens equation and use ubh instead. However, we must take 
care of the spatial dependence of Auq and Aug . We can expand them around Ubh = and write 

Au = Au^ + Jqubh (39) 
Au s = a4 0) + J s u B h (40) 

where Jo and Jg are the Jacobian matrices of Auq and Aug evaluated at Ubh = 0. 

In the same way, we can express the perturber deflection angle as in the classical Chang & Rcfsdal lens, expanding 
to first order in u 

D L s 

a Q = -pr - olq = ag(O) + Jqubh (41) 

where kq and 7q are the convergence and the shear of the perturber at the location of the black hole, respectively, 
and the orientation is chosen to diagonalize the tidal matrix Jq |2l| . 

After all these considerations on the terms of Eq. I|37|l , we can write the lens equation as 



p = £l>bh + Au^ + J u BH ] - ^A^ H (u BH ) - a Q (0) - J Q u BH + D^iAuf + J s u - Au^ - J u] (43) 



We can englobe all constant terms in /3, redefining the origin of the source plane. Rather than introducing a new 
notation, for simplicity we will continue to use /3 and drop the constant quantity A/3 from the right hand side of 
Eq. @3J, 

A/3 = ||Au^°» - a Q (0) + -^Auf ~ -a o (0). (44) 

Note that the shift in the source plane is largely dominated by olq(0), which is of order tpq, while all other shift terms 
are at most of order ipob/ D§ and we remind that we are assuming the scale of variation of the external field b to be 
much smaller than all distances among source, lens and observer. 
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Similarly, all spatial variations of the shift Auo are negligible w.r.t. the deflection angle spatial variation, so that 
they too can be neglected from the lens equation, which reduces to 



p = D L u B h ^— Aa n (ubh)-JqUbh- 



(45) 



As anticipated before, the shifts play no role in the lens equation. However, Auo intervenes in the relation between 
and itBH an d must be taken into account to describe the correct position of the images in the sky. 

Aa™ can be expressed by Eq. (JTUJl where now e = (u B h - £>l#°) / and 0° is still given by Eq. @. Writing 
the components of Ubh as D\,9^(l + e) (cos </>, sin 4>) and those of (3 as (/3i,/? 2 ), the lens equation can be written to 
zero order in e as 



fa = € 



1 — KC 
1 - Kr 



COS ( 



sin < 



(46) 
(47) 



where ej^ is given by Eq. ()12fl. 

The perturbed lens equation has an amazingly simple form, which encourages its complete investigation. Let us 
study critical curves and caustics. The determinant of the Jacobian matrix is 



J 



1 



[e- (1 - kq +7 Q cos20)e|f] 



(48) 



E 

n ■ 



Setting J — 0, the critical curves are found as 

e = (1 - k q + 7 Q cos 2<j))e. 
which, in terms of 6, correspond to 

6 = D^Au^ +e° n [l + (l- KQ + lQ cos20) e E](/ + Jo) 



COS( 

sin t 



(49) 



(50) 



where we have denoted the two-dimensional identity matrix by /. The circular symmetry of the Einstein rings is 
broken both by 7q and by Jo, which give contributions of similar relevance here. The convergence kq only acts as 
an isotropic focusing factor. Note that the perturbations to the critical curves stay small even if kq and 7q are larger 
than 1. This is quite different from the weak field Chang & Refsdal lens, where the critical curves and the caustics 
change their topology when the convergence and the shear become larger than 1. 

The perturbed critical curves are mapped in symmetric diamond-shaped caustics centered at (3 = 0, with parametric 
equations 



1 =-27Q0°: 



cos^ 
sin 3 



(51) 



The half- width of the caustics 7q 6 n decreases with the number of loops, tending to the asymptotic value 7q 8 m . In 
principle, a point-like source might be inside the first caustics but outside the caustics corresponding to large values of 
n. The shift of the caustics from the optical axis is given by Eq. (|44|) and is the same for all of them (it even coincides 
at the first order with the shift of the weak field caustic) and amounts to the zero-order term in the expansion of 
the deflection angle due to the external field qq(0). This is very different from the effect of an intrinsic angular 
momentum of the black hole 0, Ho| . where the size and the shift of the caustics increases linearly with n. So, the 
degeneracy between external shear and intrinsic spin of the lens can be broken, in principle, observing relativistic 
images of different order in n. 

The solutions of the general lens equation can be found by reducing it to a one-dimensional form. Adding the 
squares of cos</> and sin</> yields a fourth-order equation for e that can be solved by standard methods. For each 
number of loops n around the black hole, the maximum number of images is four. The expressions of the solutions in 
the general case are lengthy and we do not report them here. Basic properties of image multiplicity can be obtained 
by considering images for a source on one of the axes, for example 02 — 0. We have the following possibilities. If 
sin</> = 0, then /3i = ±#"(1 — &q — 7q — e/e^)- Two images form at <fi = and <f> = it , with 



(1 - fee. 



(52) 
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This couple of series of images always exists and matches the unperturbed solutions when the external perturbation is 
switched off. If s'mcf) ^ 0, then it must be (1 — /cq + 7q — e/e^) = 0. For a source inside a caustic, when the condition 
(3i < 27q#" is fulfilled, an additional couple of images is produced, whereas sources outside it have only two. The 
new sequences of images form at 

e = (l-fc Q+7Q ) e E (53) 
<p = ± arccos ( — - —p^. ] ■ (54) 



Note that the caustic is fully inverted w.r.t. the critical curve. This means that if the source enters from the 
top-right side, the new images form on the left-bottom side. This is different from what happens in Kerr lensing 10], 
where the caustic is only inverted on the left-right direction. 

As usual, to translate these values of e and (f> into positions in the observer sky, we must use Eq. 132fl . which gives 

= i^Aug> + + e)(J + Jo) ( ™ S f ) . (55) 



We can finally note that the total magnification of an image at position (e, (f) is given by the inverse of the Jacobian 
(|48|l . since the transformation from «bh to 9 through Eq. i|32|) introduces negligible corrections to the Jacobian 
(of order (/?q)- It is interesting to note that the two eigenvectors of the Jacobian matrix keep a quasi-radial and a 
quasi-tangential direction. We can thus distinguish a radial and a tangential magnification, which respectively read 

Mr = -e« (56) 

Mt = ~ e -(l- KQ +7QCOs20)eE- (57) 

These expressions nicely reproduce those given in Refs. 0, 0, 0, in the limit kq,7q — > 0, once we susbstitute 
the expression of (|12fl and the definition of e. Note that the radial magnification is always positive (e^ < by 
its definition), while the tangential magnification is positive if the image is outside the critical curve and negative 
otherwise. It diverges for degenerate images on the critical curve. Thus, as in all other cases in which relativistic 
images are highly magnified, they appear as elongated and very thin tangential arcs [Toj ]. 



B. Retro-lensing configuration 

The generalization of the retro-lensing equation l|15(l to the presence of an external field proceeds in the same way. 
Equation Ij32() is still valid, while Eq. (|35|1 must be revised because the source is on the same side of the lens. Starting 
from Eq. lp?4"|) . we now have xs = (DsPi, -Ds/?2, -Dls), xs = (di — #i,d2 — 02,-1) and Tg ~ —-Dls- Note that the 
z-component of the deflection vector ~ —2 is responsible of the sign flip in the z-component of the tangent vector. 

Moreover, the photon must hit the lens plane at the point — mbh in order to come back from the point Mbh- Then 
we need to integrate along the unperturbed trajectory x(°)(t) = (— uf H , — u BH , — t). Evaluating Eq. I|34|l at r = 0, 
we have 

- u BH = x(0) = D s + D LS (a ± -6)- Au s (58) 

OO OO 

Au s = 2 J dz' [\7 ±l p Q (-u¥ u ,-v% H ,z")dz", (59) 

z' 

where we have indicated by a± the projection of the deflection vector on a plane perpendicular to the z-axis. Compare 
these relations with Eqs. (|35fl and l|36|) . 

Combining Eq. I|58|l with Eq. (|32|l . we get the perturbed Schwarzschild retro- lens equation 



D s J D s 



P = _ ( i + 2^ ) - a™ + a Q + —° T ; — (60) 



Aug + Aur 



where now a BH = -AaJ| H . The minus sign is easily understood in the following way: the projected deflection vector 
increases its modulus as < q;bh < 7r/2. Then, it has modulus equal to unity at dBH = tt/2 and then starts to 
decrease as ir/2 < Abh < 7T. At q;bh = tt 5 ot^_ = 0; with a positive AaJJ H it starts to increase its modulus again but 
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its new direction is opposite to the original one. After this consideration, we can appreciate how the perturbed lens 
equation lt)U|) reduces to the unperturbed retro-lens equation (|15f) when the perturbation vanishes. 

Now let us expand the shifts Au§ and Alio for small values of u. Comparing Eq. 1)59(1 with Eq. (J2HJl, we get 

Au = Au^ + JoUbh (61) 
Au s = Au^-JoUbh- (62) 

In fact, the two shifts are calculated along very similar paths. Thus, their first order expansion in «bh has the same 
zero order term and opposite first order terms. 

Something similar happens when we calculate the deflection angle of the perturber ckq, since the incident and exit 
paths are very close each other 

oo oo oc 

a Q = 4 / V^ Q (0, 0, z)dz + 2 / d^ Q t6 B H dz - 2 / 3^ Q "bh dz, (63) 



where the symbol <9V ±<PQ/duBn indicates the Jacobian of the transverse gradient of the potential w.r.t. to «bh- 

While the zero-order contributions in the approach to the black hole and the way back sum up, the first order 
contributions in ubh cancel out. We thus conclude that, up to octopole terms, the perturbed Schwarzschild lens in 
the retro-lensing geometry is equivalent to a shifted Schwarzschild lens. Alas, the octopole contributions are of the 
same order of the contribution coming from the strong field region. Their complete investigation requires a careful 
evaluation of the perturbation of the photon trajectory in the strong field regime, unnecessary in standard lensing. 
This higher order investigation is beyond the scope of this work. So, we can conclude this section saying that, up to 
perturbations of the order ipQ(2m) /b, retro-lensing caustics are shifted but remain point-like. The shift is generally 
different from that suffered by standard lensing caustics and it is approximately given by the first term on the right 
hand side of Eq. i|tj3|) ■ As regards the images and the critical curves, everything proceeds as in the unperturbed case. 
The only difference is that 9 is given as a function of e and <fi by Eq. (|55|l . 



V. CONCLUSIONS 



Super-massive black holes are supposed to be hosted by galactic bulges, where the gravitational field of the dense 
stellar environment may perturb gravitational lensing phenomenology. While the perturbation of external fields on the 
Schwarzschild lens in the weak deflection limit is well-known, our paper has clarified the effects of a tidal perturbation 
on relativistic images generated by photons winding one or more times around a Schwarzschild black hole. 

We have firstly given a form of the unperturbed lens equation more general than that given in Ref. Q . 

By a dimensional argument, we have then demonstrated that the perturbation to the photon trajectory mainly 
arises in a weak field regime, since the perturbation in the strong regime is of higher order in the local curvature 
radius produced by the perturbing body. On the other hand, the black hole deflection in the weak field regime can be 
safely neglected, since the photon trajectory is almost face-on and the transverse gradient of the black hole potential 
is very small. 

In this way, we have practically decoupled the two deflections. This justifies the construction of the lens equation 
as a three-steps process. First we have the weak deflection by the perturber before the photon reaches the black hole, 
second we have the strong deflection by the black hole and last we have another weak deflection by the perturber. 
We have explicitly built the lens equation for the two most interesting cases: standard and retro-lensing. 

For the standard case (source behind the black hole) the perturbed lens equation reveals that the infinite sequence 
of relativistic Einstein rings near the photon sphere is slightly distorted and the point-like caustics get a finite size, 
being shifted by a constant amount. The size of the caustics decreases with the number of loops, being of order of 
the size of the shadow of the black hole multiplied by the external shear. A source inside a diamond-shaped caustic 
generates two additional relativistic images, as usual. 

In the retro-lensing geometry, the perturbation of the external field cancels out at the lowest order. In order to 
evaluate the effect, it is necessary to go to the next order in 1/6, which also involves the perturbation of the photon 
trajectory in the strong field region around the central black hole. 

The present work certainly does not exhaust all possible cases of tidal distortion, since we have restricted to weak 
perturbers in a static situation. However, the most physically relevant cases are clearly included in this class. In 
fact, the gravitational field at the Galactic center generated by the environment surrounding the supermassive black 
hole should be reasonably described in the weak field approximation. Gravitational lensing by Sgr A* could be a 
very accurate tool to track the existence of a non- vanishing shear in the field generated by the environment. In other 
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galaxies, where sometimes the bulge hosts more than one supermassive black hole, gravitational lensing may help to 
detect the mass ratios of different black holes. In any case, given the ubiquity of shear fields, the study presented here 
fills an important gap in the literature about gravitational lensing in strong gravitational fields. 
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